Role of commensurate arrangements in the optical response of metallic gratings. 
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Light localization on commensurate arrangements of metallic sub- wave length grooves is studied. 
We theoretically show that as the degree of commensuration tends to an irrational number new 
light localization states are produced. These have properties close to that reported for hot spots on 
disordered surfaces and are not permitted for simple period gratings. Some theoretical predictions 
are experimentally provided in the infra-red region by reflectivity measurements performed on two 
commensurate samples with respectively two and three slits per period. 
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Metallic systems such as random surfaces, thin films 
deposited on cold substrate or gratings, are known for 
their various and sometimes puzzling optical properties. 
In another context it is also known that quasi-crystals 
giving rise to a " forbidden" five-fold symmetry as well as 
one dimensional incommensurate structures may present 
electron localization effects [TJ [2j [3J [4]. We thus expect 
new light localization effect by merging electromagnetic 
resonances of metallic surfaces and ID incommensurate 
structures. The arrangements we study are sketched on 
fig.l. They are composed of several identical closed rect- 
angular grooves per period separated by two possible sub- 
wavelength distances L (for long distance) and S (for 
short distance). The sequence of lengths, LSSLS for in- 
stance, is chosen to be uniform and the arrangement is 
commensurate [4j|5j [6]. We can by this way study the op- 
tical properties of surfaces with complex topologies and 
interestingly we find that they present strong and sharp 
resonances associated to very local near-field intensity 
enhancements. Systems with N grooves per period, but 
separated by the same distance, have been theoretically 
studied [3 [8]. In the present context, they can be seen as 
particular cases of our commensurate structures given by 
the sequence of lengths SSS..(xN)L. Up to three cavi- 
ties per period, both approaches are equivalent. However, 
as grooves continue to be added our arrangements tend 
to a quasi-periodic (incommensurate) structure and both 
the near and far field properties become very different. 
Indeed, the commensurate arrangements present local, 
strong and wavelength dependency of the field localiza- 
tion inside the grooves. The near- field intensity distribu- 
tions are associated to sharp cavity resonances which can 
be evidenced by dips in the reflectivity. An experimen- 
tal prove of their existence is presented via reflectivity 
measurements performed in the infra-red region on grat- 
ings with respectively two and three slits per period. In 
that particular case, these measured modes also give the 
first experimental evidence of those predicted by Skigin 
et al.[9]. 
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FIG. 1: Example of a uniformly ordered commensurate struc- 
ture < 3/5 > of period D = 3S + 2L where L and S respec- 
tively stand for the long and short distances separating two 
cavities. The grooves are identical with an aperture w and 
a height h. Inset represents a SEM image of the < 2/3 > 
sample. 

The gratings of interest are characterized by a series 
{x n }i<n<Q of location of the center of the Q cavities 
contained in a period. The structures are generated fix- 
ing the origin at x\ = and using the equation[6]: 

x p+ i = x p + S x a p + L x (1 - a p ) 
with a p =[px R/Q] - [(p - 1) x R/Q] 

where the function [y] defines the integer part of R the 
number of grooves separated by the distance S and Q 
the number of grooves per period. These commensurate 
structures are called R/Q gratings, their period is D = 
Rx S + (Q — R)xL and the ratio R/Q indicates the order 
of commensurability. It is chosen to be an irreducible ra- 
tional number. For example, the successive R/Q may be 
built by using the Fibonacci series {1,2, 3, 5, 8, 13...} and 
take de values {1/2,2/3,3/5,5/8...}. For large R and Q 
the structure tends to a quasi-periodic (or incommensu- 
rate) one[6 j? D -> oo, and R/Q -> f = (y/E - l)/2 the 
inverse of the mean gold number. We have calculated 
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the optical response of these structures illuminated by 
a p-polarized plane wave (magnetic field parallel to the 
groove) following the procedure detailed in ref[ll . The 
theoretical approach is an approximated modal method 
using the surface impedance boundary conditions [10 , 
which has been several times successfully compared with 
experiments [12] when considering gold gratings with sub- 
micron geometrical parameters in the infra-red region. 
For commensurate gratings of period D the Rayleigh de- 
velopment of the field in the region (I) above the grating 
(fig.l)is: 
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where k = 2ir/X is the wave vector of the incident plane 
wave. The terms 7 m = sinO + m\/D, /3 m defined as 



= 1—7^, and R m are respectively the normalized 
wave vectors and the amplitude of the m th order of re- 
flection. In region (II) the field is determined considering 
the vertical walls as perfectly conducting and applying 
the surface impedance approximation at the bottom of 
the cavity: 
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where \l(x) = 1 for x G [—w/2,w/2] and elsewhere. A™ 
is the amplitude of the n th mode in the p th cavity with 

: (n n /k + Z)/(n n /k- 
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) 2 . The terms 



Z\ where Z = l/y/e and e the dielectric constant of the 
metal, are the reflection coefficients of the n th mode at 
the bottom of the cavity. Determining the coefficients 
A™ and R n we find the electromagnetic field in the whole 
space. Ap is obtained by solving the matrix system: 
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where we have defined the terms S^ n and F n g as: 
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FIG. 2: (a) Specular reflectivity and (b) magnetic field inten- 
sity at the interface y — at the frequency resonance located 
around 1990 cm -1 , calculated for several uniform structures 
as Rj Q tends to the inverse of the mean golden number. Sim- 
ulations were made with h — 1/xra, w — 0.3/xra, S — 0.7/xra, 
L — 1.2/jm and an incidence angle = 15°. Vertical black 
lines in (b) indicate the position of the cavities within one 
period. 



while R n is given by: 
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Figure 2a displays the evolution of the specular 
reflectivity calculated at = 15° as the order of 
commensurability tends to £. Simulations were made 
considering only the fundamental mode (n=0 in eq.l) 
and using the gold dielectric constant from ref[13 . 
For a simple-period grating, the spectrum presents 
the well-known broad dip due to the excitation of the 
Fabry-Perot like resonance inside the cavities [14 . As 
grooves are added, by considering the successive com- 
mensurate structures 1/2,2/3,3/5,5/8, supplementary 
and sharp resonances appear. For larger Q (Q > 8) 
only slight changes appear in the spectrum; the optical 
properties converge. In particular, one can see how the 
< 8/13 > and the < 55/89 > arrangements almost 
exhibit the same features. This is a direct indication 
of the self-similarity of the structures produced as 
R/Q — > £[5 1. The convergency of the optical properties 
lead to a situation where some strong and narrow 
resonances persist and are associated to strong local 



fields. These properties are also direct consequences of 
the (in) commensurate character of the structure and do 
not exist for simple or compound gratings. Fig. 2b shows 
the progressive establishment of hot spots, which are 
features experimentally observed on disordered metallic 
surfaces [T6]. as R/Q — > £. In the chosen configuration, 
the magnetic field intensity at the interface can locally 
be more than two orders of magnitude larger than the 
incident one (fig. 2b) and this corresponds to enhance- 
ments of the electric field at the interface, larger than 
three orders of magnitude. 

The occurrence of these new sharp cavity modes is 
due to the breaking of symmetry of the system. The 
pseudo-periodicity of the fields remains valid for the 
global period D such that cavities within a period do 
not see the same exciting field and this opens up for 
different configurations. A physical image can be given 
by analogy to coupled dipoles. Indeed, in the case of a 
two-groove system a complete analytical study of the 
fields expression can be performed and allows to quan- 
titatively characterize the nature of the resonances [15 : 
each resonating cavity acts as a damped oscillating 
dipole and their near-field coupling causes the splitting 
of each individual mode into a large one corresponding 
to the symmetric coupling of the cavities (—>—>) and a 
thin one, corresponding to the anti-symmetric coupling 
(—►<—). In the same manner, R/Q arrangements act as 
Q oscillating coupled dipoles, each mode corresponding 
to a specific configuration of the dipoles' orientation 
leading to more or less radiative modes and more or 
less intense near-field enhancements. Coupled modes of 
three slits with different individual frequency resonances 
were measured in the microwave regime [17]. Here, the 
Q cavities are identical but can still resonate at slightly 
different frequencies. A consequence is that the local 
intensity enhancements can be quite critically wave- 
length dependent. This property is also experimentally 
observed on disordered [16] or SERS (Surface Enhanced 
Raman Scattering) surface [18 . This point is illustrated 
in fig. 3: maps of the electric field intensity along the 
x-axis and above a unit cell of a < 8/13 > grating are 
calculated for different incident wave-number scanned 
around the strong resonance located at 1990 cm -1 . 
Within a very narrow spectral range, here 1995 to 2038 
cm -1 , the spatial localization of the electromagnetic 
field is strongly modified. For specific frequencies, one or 
few cavity (ies) act as active sites by almost individually 
concentrating very strong field intensity while the 
neighbouring ones are extinguished (fig3a and 3c). 
To validate some of these theoretical results, we evi- 
denced the existence of the sharp modes by measuring 
the specular reflectivity of gratings with one, two and 
three slits per period. The samples were prepared by 
electron-beam lithography and a double lift-off technique 
similar to ref.[19 : A Si substrate is first covered by a 
200 nm-thick gold layer and secondly by a 1 /im-thick 
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FIG. 3: Intensity maps of the electric field along the x-axis 
above one period of a < 8/13 > grating calculated at 1995 
cm -1 (a) 2013 cra _1 (b) 2031 cm _1 (c) and 2038 cra _1 (d), with 
9 — 10°. S — 0.7/xra and L — 1.2/nm h — 1/ira and w — 
0.3/ira. The "jump" of photons localization is clearly seen, 
while the incident wave number varies of only a few percent. 

SisN4 layer deposited by plasma-enhanced chemical 
vapor deposition. SisN4 walls of the width of the 
slits and separated either by the distance L=1.2 or 
1.5/im or S=0.7/im, are performed by electron-beam 
lithography and reactive ionic etching. A 650 nm-thick 
gold layer is then deposited, and the SisN4 walls are 
lifted-off by means of a HF solution. The obtained 
grooves have a trapezoidal shape (wdown = 0.36/im and 
Wtop — 0.62/im) while rectangular grooves are considered 
in the model. Small discrepancies between some nominal 
values and those used in the calculations are attributed 
to this shape difference. Grating area are 2x2 mm 2 . 
Reflectivity measurements of p-polarized impinging light 
were performed at room temperature, in dry air, in the 
1.67 to 6 jim spectral range, with a BioRad FTS60A 
DigiLab Fourier Transform photospectrometer. Spectra 
were normalized to the reflectivity of a flat gold surface. 
Figure 4 displays reflectivity measurements and cal- 
culated spectra realized at 9 = 10°. Parameters used 
in the calculation are: ^<i/i> = 0.64/xra, D = 2.2 
/im, h <1 / 2 > = 0.69/ira, S = 0.7 //m, L = 1.5 /im, 
^<2/3> = 0.6/ira, S = 0.65/im, L = 1.3/im and 
w = 0.55 jam for all samples. Figure 4a shows the 
previously measured large cavity mode noted CM as 
well as the n = ±1 branches of the surface plasmons 
noted SP[11 . These SP occurs around the same wave 
number for the two other gratings since the period of 
the three samples have close values. More interestingly, 
the splitting of the large CM into the two predicted reso- 
nances for the < 1/2 > grating and into three predicted 
resonances for the < 2/3 > grating is measured. The 
three resonances of < 2/3 > correspond to the symmetric 
coupling of the three cavities (—>—>—>), the symmetric 
coupling of the two external cavities, anti-symmetric 
to the central one whose dipole momentum is twice 



FIG. 4: Measured (black line) and calculated (light gray line) 
specular reflectivity at = 10° of the gratings < 1/1 > (a), < 
1/2 > (b) and < 2/3 > (c). Cavity resonances are noted CM 
and surface plasmon SP. The arrows indicate the direction of 
the equivalent dipole momentum at the mouth of each cavity. 
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larger (— h ►) and the ant i- symmetric coupling of the 

external cavities leading to the extinction of the field in 
the central one (— > <—). 

To highlight the symmetric or anti-symmetric nature 
of the two sharp resonances of the < 2/3 > grating, 
we have performed reflectivity measurements at various 
angles. Increasing breaks the symmetry of the system 
and favors the excitation of the anti-symmetric mode 
which vanishes at normal incidence (the symmetric 
incident field cannot couple to it). Figure 5 evidences 
the turnaround of the strength of two resonances as a 
function of the incidence angle and demonstrates that 
peak at smaller wave number is the ant i- symmetric mode 
(— » <—): it is weakly excited at small incidence angles 
and grows to become predominant around = 33°. 
Opposite behaviour is observed for the symmetric dip. 
Modifying the incident angle may thus be an easy way 
to control the near-field distribution. 
In conclusion, theoretical study of commensurate grat- 
ings tending to incommensurate structures was used to 
model and understand the optical properties of metallic 
surfaces with complex topologies. Thanks to these 
structures we can generate, in a controlled manner, 
weakly radiative modes which localize local, strong and 
wavelength dependant electromagnetic fields. These 
features are of great interest since they are commonly 
admitted to occur on disordered and/or SERS surfaces. 
Finally, experimental evidence of these modes was 
given in the infra-red region for simple commensurate 
structures. 

The authors would like to thank Nathalie Bardou for 
assistance in the fabrication process. 



FIG. 5: Experimental (black line) and calculated (gray line) 
specular reflectivity of the < 2/3 > grating performed at = 
4° (a), = 12° (b), = 24° (c) and = 33° (d). Calculations 
are made considering four modes in the cavities. 
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